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A METHOD FOR DETERMINING THE CAMBER AND TWIST OF A SURFACE TO SUPPORT A 
GIVEN DISTRIBUTION OF LIFT, WITH APPLICATIONS TO THE LOAD OVER 

A SWEPTBACK WING 

By Doris Cohen 


SUMMARY 

A graphical method is described jor finding the shape ( camber 
and twist) of an airfoil having an arbitrary distribution of lift. 
The method consists in replacing the lifting surface and its wake 
with an equivalent arrangement of vortices and in finding the 
associated vertical velocities. By a division of the vortex pattern 
into circular strips concentric about the downwash point instead 
of into the usual rectangular strips, the lifting surface is reduced 
for each downwash point to an equivalent loaded line for which 
the induced velocity is readily computed. The ratio of the ver- 
tical velocity to the stream velocity is the slope of the surface in 
the free-stream direction. 

In order to illustrate the application of the method, the shape 
of the wing consistent with the pressure distribution given by the 
two-dimensional theories has been found for three wings: an 
elliptical wing of aspect ratio 6 with 80° sweepback is treated in 
some detail; final results are presented for two unswept elliptical 
wings of aspect ratios 6 and 8, respectively, to show the degree of 
approximation involved in the use of the two-dimensional theories 
for three-dimensional flow. It is concluded that the two-dimen- 
sional theories are adequate over most of the span of a straight 
wing but that the more exact treatment is necessary in the tip 
regions or if the wing is swept. 

Application of the method to solve the reverse problem— -finding 
the lift distribution over a given surface — is also illustrated. 
The load over an uncambered, untwisted, elliptical wing with 30° 
sweepback has been calculated by successive corrections of the 
assumed vortex distribution. The result indicates a If-percent 
loss in total lift due to the introduction of sweepback with the 
greatest loss taking place at the center of the span. 

INTRODUCTION 

Because of the effect that the pressure gradients over the 
surface of a wing have on the drag, it would be of consider- 
able advantage to be able to specify the camber and the 
twist of a wing that would produce a desired distribution of 
lift. Present methods of airfoil design depend on two- 
dimensional-flow theories, which treat the spanwise and 
chordwise components of the flow independently. Although 
a theoretical treatment of three-dimensional flow is given by 
Prandtl (reference 1) and calculations have been made for 
special cases (references 2 and 3), no practical procedure for 
the arbitrary lifting surface is indicated. 

In the present paper, a method is described whereby the 
camber and the twist of a surface of arbitrary plan form may 


be determined so as to support a specified distribution of 
lift. For this method, the lifting surface and its wake are 
replaced by a distribution of vortices in a plane. The verti- 
cal velocity induced at any point on the surface by the vortex 
system defines the slope of the surface at that point. Thus, 
the problem becomes the determination of the induced veloci- 
ties. A method is presented for determining these velocities 
which, by employing chiefly graphical means, eliminates the 
difficult integrations that have limited previous work. 

The substitution of a plane vortex sheet for the lifting sur- 
face is analogous to the standard procedure of the two- 
dimensional thin-wing-section theory (see, for example, 
reference 4, p. 87), in which the flow about a thin, cambered 
section is approximated by an arrangement of vortices along 
the chord line. Inasmuch as the induced normal velocities 
are assumed to be substantially the same at the chord line 
and at- the airfoil, the ratio of these velocities to the free- 
stream velocity gives the slope of the camber line. In three- 
dimensional flow a reference plane is assumed, so situated 
that the airfoil may be considered to be a slight deviation 
from it. Upon this plane the plan form and the pressure 
distribution of the surface are projected. As in two- 
dimensional treatments, the slope is calculated in the free- 
stream direction. 

The substitution of a vortex sheet for a lifting surface is 
discussed at some length by Von Karman (reference 5, p. 15). 
In the application of the method, however, the vortices are 
generally assumed to have a rectilinear distribution. Even 
with this limitation, the evaluation of the integrals involved 
in finding the induced velocities presents considerable diffi- 
culty. (See ch. IV, sec. 15, of reference 5, where the formulas 
are developed for a rectangular wing.) The integration is 
greatly simplified by the introduction of polar coordinates 
so chosen that the elements of integration are circular strips 
concentric about the point at which the downwash is to be 
found. Application of this method is not restricted to any 
particular arrangement of vortices or to any specific form 
of the surface. 

Although the procedure described is for the determination 
of the surface that will fit a required pressure distribution, 
it may also be adapted to effect the reverse analysis, that is, 
to find the pressure distribution over an arbitrary surface. 
Such an application of the method is made to find the pres- 
sure distribution over a fiat wing of elliptical plan form 
swept back 30° at the 50-percent-chord line. 
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DETERMINATION OF THE VORTEX PATTERN FROM THE 
PRESSURE DISTRIBUTION 

The vortex pattern is obtained by integrating the chord- 
wise pressure distribution back from the leading edge at. 
several stations along the span. The circulation or vortex 
strength T will be shown to be proportional to this integral; 
the lines connecting the points where the values of the inte- 
gral are equal therefore define the vortex lines. 

In figure 1, a distribution of lift is arbitrarily specified for 
a tapered wing in straight flight. The corresponding vortex 
lines are drawn on the plan form of the wing and in the wake 
in figure 2 to show a typical vortex pattern. 




Figure 2.— Contour lines of circulation function, or vortex pattern, obtained by integrating 
the pressure distribution of figure I. 

The demonstration of the relation between the pressure 
distribution and the vortex pattern is given in the following 
paragraphs. 

In the replacement of a lifting surface by a vortex sheet, 
the assumption is made that the pressure increments due 
to the presence of the airfoil in the stream are equal and op- 
posite on the upper and lower surfaces, as would be true in 
the case of a thin plate at a small angle of attack. The 
substitution is still admissible in the calculation of lift when 
the thickness is not negligible, because the difference between 
the velocities on the upper and lower surfaces at any position 
and not the magnitude of each increment determines the 
lift at that point. 

Thus, let AV b and AV t represent for any point the local 
velocit5 7 increments on the upper and the lower surfaces due 


to the plate or to the vortices that are equivalent to it; let 
u u and u t be the magnitude of their components in the direc- 
tion of the free-stream velocity V; lot o u and t'i be the compo- 
nents normal to V in the plane of the wing; and let v.\ and 
Wi be the components normal to the wing. Then the 
pressure on the surfaces would be, if p is the air density, 



V+AV, 


pKV+UvY+vJ+Wu*] 


— 2 p( ^ 2 +2l Uu-YuY -\-vY 

If second-order effects are neglected, this expression reduces 
to 

Vu~2 p(U 2 'h2Th u ) 

Similarly, 

lh=\p{V i -\-2Vu l ) 

The resulting lift per unit area is then the difference in 
pressure, or 

p(2F)(« u — «;) (1) 

The derivation of the equivalent vortex pattern follows 
directly from equation (1). If d,T represents the element 
of circulation around a small length ds, parallel to V, over 
which the velocities u ti and u t may be considered constant, 
the following relation bolds: 

.. d s V=(u u — ui)ds (2) 

from which equation (1) may be rewritten 

dr 

Thus, the lift at every point is proportional to i or the 

cross-stream component of the vortieity. Consider now a 
narrow strip of varying width just behind the leading edge 
of the vortex sheet, such that f Lp ds is constant along the 
strip. Such a strip would represent a vortex element of 
strength f &p ds. A second vortex element could be defined 
in the same way to lie just behind the first. From equa- 
tion (3), however, 

^ J*Ap ds—T (4) 

Equation (4) defines the function F for any point on the 
lifting surface as the total circulation ahead of the point, 1 
Thus, the boundaries of the vortex element are lines of equal P. 

J The function r is numerically equal to the difference in velocity potential lietwcen the 
surfaces, for, at any point, and where and <fn are the potential functions 

over the upper and lower surfaces of the airfoil, respectively. From equation (2), 

d< j>i 

d$ ” ds ds 

or . 
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If the vortex elements are reduced in size and increased in 
number to form a continuous vortex sheet, their pattern is 
indicated, as in figure 2, by the contour lines of the function T. 
In order to satisfy the Kutta condition that there be no 
pressure difference across the trailing edge, these lines must 
leave the wing parallel to the stream velocity, as shown 
by equation (3). They then follow the streamlines down 
the wake. The drawing of these contour lines from the in- 
tegral of the pressure distribution is the first step of the 
procedure for finding the camber and the twist of the lifting 
surface. 

DETERMINATION OF THE INDUCED VERTICAL VELOCITIES 

The induced vertical velocity at. any point of the plane is 
obtained by integrating the Biot-Savart equation over the 
entire vortex pattern just described. The resulting double 
integration is reduced to a single integral by a relatively 
simple graphical procedure if a system of polar coordinates 
having its origin at the point is employed . 

Let it be required to find the downwash w at a point P of 
the vortex pattern just described. Consider a small sector 
of the plane included between two radii from P, the angle 
between the radii being d ip (fig. 3). At any distance r from P, 



the radii will cut off a small length dl of a certain vortex 
element. If thewvidth of the element in the radial direction 

ar 

is dr, the strength of the vortex element is dr. Then, 

by Biot-Savart’s rule, the downwash at P due to the small 
length of vortex will be 


dw= — -s dr dl sin 
4?r r or 


(5) 


where (3 is the angle between the vortex element and the 
radius. The length dl sin (3 is the projection of dl on the 


circumference of the circle of radius r around P, so that 
dl sin jS=r dp 

and 

du>= t"- dr dp 
Lrr or 


(6) 


The total downwash at P would be obtained by integrating 
equation (6) over O^sr— and through 360° of p. Thus, 


w= 



1 ar 

4 xr dr 


dip dr 


(7) 


In order to evaluate the double integral of equation (7), the 
relation (Liebniz’s rule) 


X 


dr d 
d ¥ a * = ct? 



( 8 ) 


is used to reduce the first integration to a graphical procedure. 

The first step in evaluating w is to draw on the vortex pat- 
tern circles of various radii about P. Around any one circle 
(see fig. 4) the function r will take on values indicated by the 



intersections of the circle with the contour lines of T. If 
these values of T are plotted against the angle \p — measured, 
let us say, from the free-stream direction — graphical inte- 
gration will give J* r dip for each circle. It is somewhat 

more convenient to plot T against pj2ir. Then the integral 
will be the average value of T around the circle, designated 
in the usual way by 

(9) 


Then r is a function of r. (See fig. 5.) 

From equation (8), 

f 2 ' dr , 9 dr 

Jo 5 r d *- 2 *di 

(10) 

so that 


r° i df , 

w= — fr dr 

Jo 2r dr 

(11) 
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Fiocre 5.— Typical curve forf (r). (Curve for point P of fig. 4.) 


For the evaluation of equation (11), first plot r against 
r, as in figure 5. The curve will approach an asymptote, 
which is readily found from the expression for the load 

dv 

curve across the wake. Thus, when r is large, d\p=- 7 and 
equation (9) reduces to 

? =2UJI r <»> * < 12 > 

* 

where y is the spanwise dimension and b / 2 is the semispan 
of the wing. If the area under the load curve across the 
wake is A, then for very large values of r 



The plot of T(r) should be carried out to a value of r such 
that the curve approaches this asymptote within the 
accuracy of the work. 

The load curve of figure 5 is a typical one. The following 
method has been found particularly suited to the determina- 
tion of the downwash at the origin of such a curve. The 
first section of the curve (designated by I), starting with 
zero slope, is approximated up to the inflection point, or to 
a point r 0 somewhat ahead of it, by an expression of the 
form r=a 0 — a 2 r 2 — a 4 r 4 . Additional terms might be used 
but are generally not necessary. The downwash due to a 
curve of the general form 

f=a 0 —a a r n (n> 1 ) 


at r =0 is given by 


a n n 
2 71—1 


f rp 71 1 m 7t — IN 
V B >A 


where r A and r B are the end values of r for the interval over 
which the curve extends. The downwash contributed by the 
first section of the load r=a 0 — a 2 r 2 — a 4 r 4 is therefore 


w(I) = a 2 r 0 +g a 4 r 0 3 


The part of the curve immediately following r = r 0 lias a 
critical effect on the value of the downwash, at the same time 
being usually too irregular to be approximated for any 
distance by a simple algebraic expression. It is therefore 
advisable to proceed as in numerical integration, dividing 
the curve into a finite number of small sections and consider- 
ing each section of the curve to have a simple mathematical 
expression. Because the effectiveness of the variation of 
load depends on its closeness to the downwash point, the 
intervals are taken in geometric rather than arithmetic 
progression. Thus, the abscissas are r 0 , At 0 , A'%, . . . etc., 
where the ratio k is a number, usually between 1 and 2 , 
determined by the size of the intervals required for accurate 
representation of the curve immediately following r 0 . The 
usual procedure now would be to assume the curve to be a 
straight line over each small interval; but- when the curvature 
is largely in one direction, as it is in these curves, this assump- 
tion introduces a small but cumulative error, which may 
amount to 10 percent or more in the total. The following 
method, which fits the curve with a succession of parabolas, 
is found to give very good accuracy with no increase in 
computation. The method is best presented in tabular 
form: 


Computation q£ v> 


rjtt rjnfi 

—k in r& ^krin 


A I'ln AHhM 


ATi* 

rj» 


AHmi 

rim 


N- 1 


To 


n 


r»-ro 


fi-Fi 


rt-n 


r* 


n 


n 


r 2 .V -1 


r*y 


ri-n 


I'i-v-t — IYv-i 


ATb 
r 0 


AT* 

rt 


An 

Ti 


Arjy-t 


Kt—lc (t; 


77 logfe 
Kl “ l 


Ah 

r» 


ATi 

ri 


ah 

ri 


ATty-i 

r*v-i 


Totals: V — Hilt 1 

^-L Tu Tin 

n — Q 0 




k (i 


(14) 
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In the following table are given the constants IQ, K u and 
for several convenient values of k. Since the choice 


k(k— 1) 

of k is not critical, the values included should serve with- 
out interpolation. 


k 

K, 

-El 

1 

i(t-I) 

1.02 

0.01000 

0.01000 

49. 02 ' 

1.05 

.02582 

.02420 

19.04 

1. 07 

. 03567 

.03340 

13.35 

1. 10 

^05325 

.04690 

9.091 

1.15 

.08222 

-0GS27 

5.797 

1.20 

.1127 

.08S40 

4.167 

1.25 

. 1446 

.1074 

3.200 

1.30 

.1780 

. 1255 

2.564 

1.50 

.3246 

.1891 

L 333 


The value of the ratio k that will give sufficiently small 
increments of r where the slope is large will probably be 
found to be smaller than necessary after the curve has become 
more regular and the distance from the origin larger. The 
computations may be interrupted here and the downwash 
10 ( 11 ) contributed by the second section may be calculated. 
Then, using the last abscissa r 2y as a new starting point and 
a larger value of k, compute in the same way the downwash 
ifl(III) due to the remainder of the curve to a point where 
the difference between the curve and its asymptote is 
negligible. 

For the portion of the curve extending to infinity (section 
IV of fig. 5), the previously determined asymptote is used 
and the downwash found analytically. Thus, from a large 
value R of r out to infinity, 


w(IV) = 



i dr 
r dr 


dr 


from equation (10) ; and if dT/dr is found from equation (13), 

w(IY)= db ( 15 ) 

The downwash at P is then the sum 


10(1) 4- 10(11) 4- wflll) 4- 10(IV) 


It is interesting to note that the three-dimensional problem 
has been reduced to one of two-dimensional flow, as may be 
seen by replacing r with its original expression (equation (9)) 
in equation (11), which may then be written 



(16) 


Equation (16) is recognized as the usual formula for the 
induced normal velocity, with the load expressed in the 
form of a definite integral, and suggests that the first integra- 
tion (except that the factor 1/2 ?r was introduced) was equiv- 
alent to concentrating all the vorticity around each circle 
at a single point at the distance r along a line of infinite 
extension from P. The loaded line of figure 5 may be con- 


sidered, except for the factor 1/2 x, to be the equivalent of the 
original lifting surface in its effect at P. 

EXAMPLES 

The method will be applied to check the elliptical dis- 
tribution of lift conventionally assumed for an uncambered 
elliptical wing. This distribution, arrived at by combining 
the two-dimensional theories, does not take account of 
sweepback or stagger of the lifting elements. In the ap- 
plication of the present method, three cases of elliptical 
chord distribution will be investigated, two with a straight 
50-percent-chord line, and one with the 50-percent-chord 
line swept back about 30°, the sections remaining parallel 
to the plane of sy mm etry. The calculation of the vertical 
velocities over the sweptback wing will be carried out herein 
in some detail in order to illustrate the method. The down- 
wash will be found at several points along the three-quarter- 
chord line. The downwash at the three-quarter-chord 
point of a section is of particular interest because, from the 
thin-wing-section theory (reference 6, p. 82), the effective 
angle of attack of the section is given by the slope of the 
camber line in the neighborhood of that point, if the camber 
is approximately circular. 

SWEPTBACK WING 

For the chord wise lift distribution, the two-dimensional 
flow around a flat plate as given by the thin-wing-section 
theory is assumed. This lift distribution and the circulation 
function obtained by integrating it are shown in figure 6. In 



Figure 6.— Cfcordwise pressure distribution and circulation function assumed for figure 7 
from the two-dimensional-how theory for a fiat plate. 




548 


REPORT NO. 82 6 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


tliis case the mathematical expression for the pressure differ- 
ence is known and r can be found analytically. For con- 
venience, the units have been so chosen that the maximum 
value of P is 1.0. If, further, all lengths are expressed in 
terms of the semispan (measured perpendicular to the plane 
of symmetry) as the unit of length, the total circulation 
(r at the trailing edge) is simply where y is the distance 

spanwise from the center line. The resulting contour lines 
of F appear, for the sweptbaek wing, as in figure 7. The 
points at which the downwash will be found are also shown. 



Figure 7.— Contour linos of elllptically distributed load over sweptbaek wing and wake. 

Aspect ratio, 6. 


In order to find the downwash at point B, circles spaced 
as shown in figure 8 are drawn about this point, and the 
values of F indicated by the contour lines intercepted by 
each circle are plotted against the angular location ip of the 
points of intersection. Of the curves for Tty) corresponding to 
the circles of figure 8, five typical ones are shown in figure 9. 
The curves designated for r=0.067 and r=0.20 are 
characteristic of circles close to the downwash point; the 

curve for r=0.216 includes a point ( =0.883') at which the 



Figure 8. —Circles of integration drawn on vortex pattern. The circles corresponding to 
the curves of figure 9 are identified. 


circle is tangent to the leading edge of the wing; the circle 
of radius 0.733 lies partly ahead of the wing, where r 0, 
partly on the wing, and partly in the wake; life circle 
of radius 1.60 traverses the vortex pattern across the wake 
only. 

The results of integrating these curves are plotted against 
r in figure 10. Since the circulation across the wake is tyl — if, 
the integral A is tt/ 2 and the asymptote for the curve of F(r) 

is, from equation (13), r=~ The downwash is now calcu- 
lated as follows: 

It is found that the polynomial 

F=0.8335-2.70r 2 -C.57r 4 

fits the curve of F through r=0.20. Then the downwash con- 
tributed by the section from 0 to 0.2 is, from equation (14) 

w(I) = 0.575 

Section II, with k small, is taken to include both inflection 
points (there is a third inflection point at 1.78, but its effect 
is negligible) . The following table shows the downwash com- 
putations for this section and for the succeeding one: 


r 

f 

AT 

A fu 
TU 

a7».„ 

TU 

Section II; k— 1.07 

0. 200 
.211 
.229 
.245 
.262 
.280 
.300 
.321 
.344 
.368 
.393 

0.715 

.690 

.657 

.636 

.620 

.606 

.593 

.552 

.569 

.560 

.550 

-0. 025 
-.033 
-.021 
-.016 
-.014 
-.013 
011 
-.013 
-.009 

010 

-0.125 
-.002 
-.053 
-. 033 
-.020 

-0. 165 
070 
050 
-.013 
-.029 

Totals: 


stem “0.320 

-.329 

-.357 

Section III; £=1.20 

0. 393 
.472 
.566 
.679 

. 815_ 

.978 
1.173 
1. 103 
1.790 
2. 028 
2.433 

a 550 
.524 
.497 
.453 
.379 
.317 
.268 
.219 
.153 
.131 
.108 

-0.026 

027 

-.044 
-. 074 
-.062 

Q49 

-.049 
-.066 
-.022 
- 023 

-0.000 
— . 07S 
-.070 
-.012 
-.012 

-o.oco 

-.131 

-.000 

-.050 

-.013 

Totals: 


W(III) =0.250 

— . 274 

-.329 


Section HI extends to r=2.433, where r=0.1Q8 and 
j-=0.1037 The downwash induced by the asymptote from 


r=2.433 out to infinity is only 


1 

16(2.433) 2 ' 


or 0.011, about 1 


percent of the total due to the curve of F, and the 5-percent 
error in the ordinate may therefore be neglected. The total 
downwash at B is then 


w(I)+7c(II)+w(III) + ro(IV) = 1.15G 
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r=0.067 


Figure 9. — Typical curves for the variation of r around a circle. The curves shown correspond to the circles identified in figure S 


T t sem/spans 

Figure 10. — Curve of f (r) for point B of figure 7. 


The downwash at points A, C, and D is found in the 
same way. Figure 11 shows the curves of r for these 

points/ The asymptote r=-^ is, of course, common to ail 

the curves for this wing. 


The downwash is plotted against the spanwise location of 
the points in figure 12. The quantities T maz and b/2, hereto- 
fore assumed to be unity, are included to make the result 
non dim ension ah This curve of the downwash or vertical 
velocity at the three-quarter-chord line, since the slope of the 


843107 — 50 36 
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surface is w/V, is a measure of the amount of twist the wing 
must have to sustain the assumed distribution of lift. If the 
wing were actually flat, which was the premise in deriving 
the distribution from the two-dimensional theories, w would 
be equal to 1.0 all over the surface, so that the devia- 
tion of the curve from the line w=l indicates the amount 
by which the two-dimensional theories are in error when 
applied to three-dimensional flow. The discontinuity in 
vorticity at the center line gives rise to a discontinuity in the 
downwash, which goes to infinity everywhere along the 
center section. This result indicates that the assumed condi- 
tion, in which the vortex lines bend to form an angle, cannot 
exist in practice. 




liNSWEPT WINGS— ASPECT-RATIO EFFECT 

The corresponding curve for the straight elliptical wing of 
aspect ratio 0 is shown in figure 13. In this case, the down- 
wash was computed also at points along the quarter-chord 
line. The slope of the surface was found to be less than the 
slope at the corresponding three-quarter-chord points. The 
positive camber thus indicated is very small at the center but 
increases sharply near the tips. 

The corresponding curves for the wing of aspect ratio 3 
are shown in figure 14. Although the deviation from the 
lifting-line-theory value is greater, on a percentage basis, 
than in the case of aspect ratio 6, the agreement with lifting- 
line theory is still rather good. It may be concluded that, 
except near the tips, the two-dimensional theories may bo 
considered generally adequate for straight wings. 



Figure 12,— Vertical velocity w at the three-quarter-chord line of the swept back wing with 
distribution, of load calculated by two-dimensional-flow theories. 


Figure 13.— Vertical velocity w at points on straight elliptical wing, aspect ratio of 0, with 
distribution of load calculated by two-dimensional-flow theories. 
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Figcee 14. — Vei deal velocity a at points on straight elliptical wing, aspect ratio of 3, with 
distribution of load calculated by two dimensional-flow theories. 

AN APPLICATION TO THE REVERSE PROBLEM: THE LOAD 
OVER A FLAT, SWEPTBACK WING 

The results of the preceding calculations for the swept back 
wing (fig. 12) indicate that for a flat surface the vortex lines 
over such a wing would actually be rounded off at the center 
and the load in that region thereby reduced below that of the 
adjoining sections. It is also apparent that the downwash 
induced over the tip regions would have to be increased; 
this would be accomplished by increasing the density of the 
vortex lines (or lift) toward the tips at the expense of the lift 
farther inboard. 

With these requirements in mind, a second approximation 
to the lift distribution was assumed; and the downwash at 
points along the quarter-chord and three-quarter-chord line 
was found. This approximation proved to be in error in 
the opposite direction; that is, the downwash now increased 
from center to tip. A slight modification of this second 
distribution made the resulting downwash variation linear 
along the span with, however, a smaller slope than given by 
the two-dimensional distribution. Values of w for the quarter- 
chord points fell along a line parallel to that for the 
three-quarter chord and approximately 8 percent below it, 
indicating a small amount of camber. 

It was assumed that interpolation between the second 
distribution, as modified, and the first (two-dimensional) 
approximation, at the same angle of attack, would be a 
fairly accurate solution to the problem. The curves pre- 
sented are the result of this interpolation. 

Figure 15 shows the final configuration of vortices for the 
flat sweptbaek wing. The concentration of lift in any 
region is proportional to the density of the lines. The entire 
pattern is independent of angle of attack, except as the basic 
theory breaks down at large angles of attack. 

In figure 16 is shown the span loading. The calculated 
load is compared with the elliptic load, which has been found 
to be fairly accurate for the straight wing at an aspect ratio 


of 6. The induced camber noted in figure 13 for this case 
is approximately the same, on the average, as remained in 
the sweptbaek wing after the last adjustment, so that the 
wings may be considered comparable. 

When the two wings are at the same angle of attack, 
measured by the slope of the surface at the three-quarter- 
chord line, the area under the lift curve for the sweptbaek 
wing is 86 percent of that under the ellipse, indicating a loss 
due to sweepback of 14 percent of the total lift. The effect 
of sweepback on the spanwise variation of lift is indicated 
by the curve for the straight elliptic wing supporting the 
same total lift. The effects of sweepback are seen in the 
higher concentration of load at the wing tips and the depres- 
sion in the load curve near the center. 

USE OF THE METHOD AT WING TIPS 

It is expected that the method outlined herein will be 
especially useful in investigating the nature of the flow near 
the wing tip, where two-dimensional approximations no 
longer can be applied. The more accurate three-dimensional 
treatments available are also unsuitable for this purpose 
because the calculations fail to converge at the tips. The 



r=0.893 


Figcee 15 —Distribution of vorticity or circulation over fiat elliptical wing, aspect ratio of 6, 

swept back 30°. 



Figuee 16. — Span-load corves for a wing with elliptical chord distribution, showing tbe 
effects of 30° sweepback. 
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method of calculation described in this paper presents no 
particular difficulty in these regions. As a test of the ac- 
curacy obtainable, the induced downwash was computed at 
a point near the edge of a circular plate in nonlifting poten- 
tial and flow was found to check almost exactly with the 
known solution. 

Nevertheless, results obtained for the tips have qualitative 
rather than quantitative value (except for low angles of 
attack) . The validity of the theory is actually limited by the 
existence of strong tip vortices, which may cause the vortex 
sheet to curl up out of the plane in which it is assumed to lie. 
On the other hand, the high concentration of vortieity as- 
sociated with this effect adds appreciably to the drag of the 
wing, so that even a general indication of such a concentra- 
tion of load at the tips is of value. It should be possible to 
design, by the use of the present method, a wing tip that 
would avoid this effect by providing a fairly gradual tapering 
off of the load span wise. At the same time, a favorable 


chordwise gradient could be specified. Thus, it appears likely 
that an optimum tip for low drag could be deduced. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 16, 1942. 
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